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The available data on neutron scattering were analysed to constrain a hypothetical new short- 
range interaction. We show that these constraints are several orders of magnitude better than those 
usually cited in the range between 1 pm and 5 nm. This distance range occupies an intermediate 
space between collider searches for strongly coupled heavy bosons and searches for new weak macro- 
scopic forces. We emphasise the reliability of the neutron contraints in so far as they provide several 
independent strategies. We have identified the most promising way to improve them. 



PACS numbers: 03.75.Dg, 03. 75. Be, 13.40.Gp 



I. INTRODUCTION 

The existence of other forces in nature, mediated by 
new bosons, has been extensively discussed in the litera- 
ture, given their possibility in many of the extensions to 
the standard model of particle physics New bosons 
for example are predicted by most of the Grand Uni- 
fied Theories embedding the standard model, with a cou- 
pling constant of ~ 10~^. These strongly coupled bosons 
would have to be heavier than « 1 TeV if they were not 
to conflict with present observations; heavier bosons will 
be searched for at the Large Hadron Collider. Lighter 
bosons could however have remained unnoticed, provided 
they interact weakly with matter. Such bosons would 
mediate a finite range force between two fermions: 
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where g is the coupling constant, Qi and Q2 the charges 
of the fermions under the new interaction, and the range 
of this Yukawa-like potential A = is inversely propor- 
tional to the boson mass M. In the following we consider 
the interactions of neutrons with nuclei of atomic number 
A: the charge of the atom under the new interaction is 
equal Qi — A; the neutron charge is equal unity Q2 — 1. 
A new boson could even be massless, as has been sug- 
gested by Lee and Yang well before the birth of the 
standard model, to explain the conservation of the baryon 
number. This additional massless boson would mediate 
a new infinite-range force, and could be seen in searches 
for violation of the equivalence principle. The presence of 
very light bosons (M ^ 1 eV) would be shown by devia- 
tions from the gravitational inverse square law. Gravity 
has been probed down to distances of 0.1 mm Q; new 
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bosons lighter than 2 x 10^"^ eV must thus have a cou- 
pling constant lower than the gravity strength between 
nucleons, < 10^'^^. 

Theories with extra large spatial dimensions [1, H, 0, 
B 3 provide strong motivation to search for such forces. 
If a boson is allowed to travel in large extra-dimensions, 
with a strong coupling constant in the bulk, it behaves 
in our 4D world as a very weakly coupled new boson, 
the coupling being diluted in the extra-dimensions. The 
light dark matter hypothesis also argues in favour of the 
existence of new short range interactions . 

While measurements of the Casimir or Van der Waals 
forces (for a review, see e.g. [ll|) give the best con- 
straints in the nanometer range (10 nm < A < 1 /im), 
and antiprotonic atoms constrain the domain below 1 
pm [ll] , it has been suggested that experiments with 
neutrons could be competitive in the intermediate range 
O Q d, d, [S mi • Neutrons could also probespin- 
dependent interactions in a wider distance range [l9| , or 
spin-independent interactions in the range of several mi- 
crometers [il,!!!. 



In this article we give the quantitative constraints on 
the parameters of the additional interaction, A and g, 
from the existing data on neutrons scattering at nuclei. 
In section|lI]we analyse the influence of a new short-range 
interaction on the scattering of neutrons at nuclei. In sec- 
tion [nil we use the fact that the nuclear radius, as well 
as the scattering lengths, are expected to be proportional 
to A^^^, where A is the number of nucleons, whilst the 
contribution of an additional interaction would result in 
an additional linear term in the mass dependence of the 
scattering lengths. In section IIVI and |Vl we use the dif- 
ferent sensitivities of different types of neutron scattering 
experiments to extra interactions (forward and backward 
scattering) in order to constrain them. In section IVll we 
propose a way to improve these constraints p^ . 
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II. SLOW NEUTRON / NUCLEI INTERACTION 
WITH EXTRA-SHORT-RANGE INTERACTIONS 

The scattering of slow neutrons on atoms is described 
by the scattering ampUtude /(q); this can be represented 
by a sum of a few terms [21| : 



/(q) = /nuci(q) + /ne(q) + /y(q) 



(2) 



The first and the most important term represents the 
scattering due to the nuclear neutron- nucleus interaction. 
At low energies discussed in this article, it is isotropic and 
energy-independant, because the nuclear radius is much 
smaller than the wavclcgth of slow neutrons: 



/„„c:(q) = -b. 



(3) 



The coherent scattering lenght b is the fundamental pa- 
rameter describing the interaction of slow neutrons with 
a nucleus . 

The second term is the amplitude of so-called electron- 
neutron scattering due to the interaction of the neutron 
charge distribution with the nucleus charge and the elec- 
tron cloud. This amplitude can be written as 



fneici)=-bne{Z-f{Z,Cl)), 



(4) 



where /(Z, q) is the atomic form- factor measured in 
the X-rays experiments and 6„e is a constant called the 
electron-neutron scattering length, which is directly re- 
lated to the neutron charge radius 21] and to the neutron 
electromagnetic form- factor Ge{(1^) by 
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m and me being the neutron and electron masses, ao the 
Bohr radius. This contribution to the total scattering 
amplitude is as small as a per cent for heavy nuclei. 

In the presence of a new interaction (HI), the scatter- 
ing for a center of mass momentum Hk due to the extra 
interaction, within the Born approximation, is given by 



fv(0) = -A^hc '-^ 
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where q — 2fcsin(6'/2), is the scattering angle. 

Any other possible contributions to the scattering am- 
plitude /(q), due to non zero nuclear radius, nucleon 
polarizability, etc. are very small in the energy range 
discussed here [HI] and have therefore been omitted in 
©■ 

The nuclear scattering lengths are measured for almost 
all stable nuclei, using a variety of methods. A review of 
the different methods and a complete table of the mea- 
sured scattering lengths can be found in [2^. We can 
distinguish two classes of method, with different sensi- 
tivities to a new interaction. 

The first class - including the interference method, 
the total reflection method, the gravity refractometer 



method - measures the forward scattering amplitude 
/(q — 0). These methods actually measure the mean 
optical potential of a given material, called the Fermi 
potential, due to the coherent scattering of neutrons at 
many nuclei. The Fermi potential is related to the for- 
ward scattering amplitude. 

In the presence of the new force, the measured scatter- 
ing lenght can be separated into a nuclear and an addi- 
tional terml29l: 
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The second class of method - including the Bragg diffrac- 
tion method and the transmission method - uses non- 
zero transferred momentum. In the Bragg diffraction 
method, the scattering amplitude for a momentum trans- 
fer of (7bd = 10 nm~^ is measured. One actually extracts, 
besides the nuclear term, an extra contribution according 

to m 



6bd = + A 
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In the case of the transmission method, the total cross- 
section is measured. Generally, neutrons with energies 
of about 1 eV are used; they are much faster than slow 
neutrons, and no coherent scattering can be observed. An 
additional interaction would manifest itself by an energy 
dependance of the extracted scattering length 
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Finally, we should also mention the very popular Chris- 
tiansen filter technique; this measures relative scattering 
lengths, so we do not consider this data. 



HI. RANDOM POTENTIAL NUCLEAR MODEL 

A simple and robust limit on the additional Yukawa 
forces can be easily obtained by neglecting the small term 
due to the neutron-electron scattering and by studing the 
general ^-dependence of the scattering amplitude. In the 
domain of A < l/q-Bu, the optical and Bragg diffraction 
methods are sensitive to the same amplitude 
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as clear from ([7]) and ([5]). The presence of additional 
forces would be apparent from the linear increase of the 
measured scattering length as a function of A. This de- 
pendence is not observed in the experimental data pre- 
sented in fig. [TJ Strong variations with A of the measured 
scattering length are mostly due to the shell effects in 
neutron-nucleus nuclear amplitude b. 

To establish a quantitative upper limit on this addi- 
tional term, the A dependence of the nuclear scattering 
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FIG. 1: Measured scattering lengths as a function of nucleus 
atomic number. 
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FIG. 2: This histogram shows the distribution of measured 
scattering lengths normalized to the radius of the nuclei. The 
curve corresponds to the random potential model. 



length b = b{A) needs to be taken into account. The 
complete ah initio calculation of these shell effects for b 
is particularly complicated and has never been reported 
in the literature. Fortunately, there exists a very sim- 
ple and elegant semi-phenomenological approach that de- 
scribes these variations It assumes that a nucleus 
can be presented as an attractive " square well" potential, 
with radius RA^^^ and depth Vq for slow neutrons. The 
scattering length would then be equal to 



biA) = RA^^^ I 1 



tan(X) 



(11) 



where X — 



1/3 



VSmVo is supposed to be a random 
variable distributed uniformly over the range [7r/2, 57r/2] ; 
the lower value corresponds to the appearance of a bound 
state and the upper limit is set sufficiently large not to 
influence the results of the present analysis; more details 
can be found i n [25| . 

Since 1970 ([25j), the number of nuclei with measured 
scattering lengths has doubled and the level of precision 
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FIG. 3; Maximum likelihood analysis of the two parameters 
R and ftExtra of (fT2)) . 



has been improved. We include here the new interaction 
analysis. 

Without any additional interaction, this analysis pro- 
vides the results presented in fig. [5) One may conclude 
that the model describes well the distribution of all ex- 
perimental data well; the value of the only free parameter 
in this model is estimated to be i? = 1.44 ± 0.05 fm at 
the 68 % C.L. The likelihood function at its maximum 
satisfies ln(L) = —254 for 216 degrees of freedom. 

With a short-range new interaction included in the 
analysis we have to consider the random variable 



RA^'^ 1 



tan(X) 
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(12) 



where the effect of the extra interaction is the slope 
&Extra — ff^A^ of the linear term. The estimation 
of the two parameters R and &Extra from the experimen- 
tal data, again using the maximum likelihood method, is 
presented in fig. [31 The linear term is compatible with 
zero, as expected. We thus obtain a quantitative con- 
straint for the coupling g{X): 



g^}? < 0.016 fm^ at 95% C.L. 



(13) 



This result is presented in fig. [7] for the distance range of 
interest, 10"^^ - 10"^° m. 
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IV. CONSTRAINT FROM COMPARAISON OF 
FORWARD AND BACKWARD SCATTERING OF 
NEUTRONS 

Another way to constrain on aditional Yukawa forces 
consists in comparing the scattering lengths measured by 
different methods. 

As explained above, the scattering lengths measured 
using the Bragg diffraction method 6bd and the interfer- 
ence method 6„pt do not show the same sensitivity to a 
new short-range interaction. According to ^ and ([8]), 
the ratio of the two values should deviate from unity in 
the presence of an additional interaction 



b 2Trhc 
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We found a set of 13 nuclei for which both measurements 
exist. 

The different measurements quoted in the litterature 
are sometimes incompatible, even for the same measure- 
ment method, because of additional systematics not in- 
cluded in the quoted errors. To compensate for this 
we estimated a methodological error a for each method. 
By selecting all the nuclei for which several measure- 
ments with a given method are available, we obtained 
the values hA,i ± A&yi.i, where A is the nucleus index and 
1 ^ i ^ riA lists the experiments using the given method. 
The methodological error is then calculated so that: 



E 



bA) 



'A,i 



A 



That is, we force the deviation from the weighted av- 
erage bA to be equal to the number of degrees of freedom 
^^nA. For the Bragg diffraction method, we have 72 
degrees of freedom, and a methodological error of 0.16 
fm; this can be compared to the average of the publised 
errors of 0.23 fm. For the interference method, we have 
22 degrees of freedom; we found a methodological error of 
0.05 fm, while the average of published errors is 0.05 fm. 
Table U summarizes. 

With the enlarged errors, a fit of eq. (|14p using the 
13 set of data for bi and 6bd, we obtain the constraint: 



{qXf 



< 0.0013 fm^ at 95% C.L. (16) 



1 + (gA)2 

corresponding to the bold limit in fig. [71 



V. ELECTROMAGNETIC EFFECTS 

Up to now, the amplitude due to a new additional in- 
teraction /v(q) has been compared to the nuclear one 
/„uci(q) (see ([2])). One could improve the limit obtained 
by comparing the amplitude due to a new additional in- 
teraction to a smaller amplitude due to an electromag- 
netic interaction (/ne(q))- This idea was first proposed 
in ref. [11. 



TABLE I: Existing data on neutron scattering lengths. 
Method Bragg Interference 

Diffraction 

Number of measurements 141 41 

Number of measured nuclei 98 28 
Mean relative accuracy 3 % 0.4 % 

Methodological error [fm] 0.16 0.05 



To achieve this, three independent measurements are 
required (roughly speaking, to determine independently 
the three contributions to (l2|)). 

One could repeat the previous analysis using measure- 
ments of the total cross-section instead of the Bragg 
diffraction. As can be seen from Q, the scattering 
length extracted at energies of ~ 1 eV {k — 200 nm~^ = 
1/5 pm) contains an extra force contribution that is dif- 
ferent from that in optical methods. In particular, if 
the range of a new interaction is larger than 1 pm, the 
scattering length extracted from the total cross-section 
at 1 eV can be considered free of any extra contribu- 
tion. The sensitivity is so high, however, that the resid- 
ual electromagnetic effects due to the neutron square 
charge radius can mimick an extra-force contribution in 
the quantity &(1 eV) — feopt, as this contribution is energy- 
dependent and proportional to the charge number of the 
atoms. This problem is known as the determination of 
the neutron-electron scattering length bne- The extracted 
(15) difference 6(1 eV) — 6opt therefore contains the two con- 



tributions: 



6(1 eV)-6(0) ^ Zb„ 



~A——g^\^ 
Zirhc 
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(17) 

Unfortunately, there is very clear disagreement be- 



tween the two groups of values for 6™p — ^ 



_ fc(l cV)-6(0) 

known as the Garching-Argonne and Dubna values [2^ 

bTf ^ (-1.31 ±0.03) X 10"^ fm [Gartching-Argonne] 
b^^P = (-1.59 ±0.04) X 10"^ fm [Dubna] (18) 

The discrepancy is much greater than the quoted uncer- 
tainties of the experiments and there evidently an unac- 
counted for systematic error in at least one of the exper- 
iments. 

In order to overcome this difficulty we could determine 
bne from the experimental data on the neutron form fac- 
tor The simplest way to do this consists in using a 
commonly accepted general parametrization of the neu- 
tron form factor [271: 



1 + 6t 



Gd, 



(19) 



where /i„ = — 1.91/iB is the neutron anomalous magnetic 
moment, r q^/4m'^ and 



Gi5(q') = 



1 



(l + q2/0.71 (GeV/c)2)" 



(20) 
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NEUTRON ELECTRIC FORM FACTOR 
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FIG. 4: The neutron form factor Gsiqi^) as a function of the 
momentum transferred q^. The experimental data are taken 
from j28(l ; the solid curve is a two parameter fit using formula 



FIG. 5: The asymmetry of neutron scattering at an atom 
in the center-of-mass reference system is shown as a function 
of various characteristic ranges of the additional interaction 
0, 1; 1; 10 nm. 



is so-called dipole form factor ; a and b being fitting pa- 
rameters. 

A fit of an existing set of the neutron form factor ex- 
perimental data [1^ yields the following values for the 
parameters: 

a = (0.77 ±0.06) 
b = (2.18 ±0.58) 

with xVNDF = 15.3/27. The results of the fit are pre- 
sented in fig. m 

Let us note that the momentum transferred in these 
experiments is very large and a contribution from the 
term /v(q) is negligible. The 6„e determined in this way 
is 



bne = (-1.13 ± 0.08) X 10"^ fm. 



(21) 



It does not agree with the value of bne , obtained in Dubna 
experiments. 

Our principal conclusion consists in the observa- 
tion of (underestimated) systematical uncertainties in 
the presented experiments. Therefore a single experi- 
ment/method can not be used for any reliable constraint. 
A conservative estimate of the precision of the bne value 
could be obtained from analysing the discrepancies in 
the results obtained by different methods; it is equal to 
Abne ^ 6 X 10^* fm. The corresponding contraint at the 
2(7 level 
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(22) 



is represented by the dot-dashed line in fig. [71 



VI. ASYMMETRY OF SCATTERING 

As is clear from fig. [3 the best constraint was obtained 
from the analysis of the energy dependence of the neutron 
scattering lengths in the bne measurements. However, 
the precision here is limited by the correction for the bne 
value itself. An obvious proposal for improving this con- 
straint would be to set up experimental conditions free 
of the bne contribution. This is indeed possible, because 
neutron-electron scattering is essential for fast neutrons 
only, and is absent for slow neutrons. 

We propose improving the experiment [2^ and measur- 
ing the forward-backward asymmetry of the scattering of 
neutrons at atoms of noble gases, in the following way: 
the initial velocity of the neutrons should correspond to 
the range of very cold neutrons (VCN); the double differ- 
ential measurement of neutron velocity before/after scat- 
tering should be used to calculate the transferred momen- 
tum for every collision. 

Fig. [5] shows the asymmetry of neutron scattering at an 
atom in the center-of-mass reference system, as a func- 
tion of various characteristic ranges of the additional in- 
teraction. It is clear that neutrons with velocities of 
1 — 1000 m/s need to be used. This means that the 
thermal motion of atoms cannot be neglected, and that 
the neutron velocity in the laboratory fixed reference sys- 
tem is not equal to that in the center-of-mass reference 
system. Even totally isotropic scattering in the center- 
of-mass reference system would be highly anisotropic in 
the laboratory system. Nevertheless, the kinematics of 
the scattering process can be reconstructed precisely if 
both the initial and the final neutron velocity in the lab- 
oratory reference system are measured. This can easily 
be done with reasonable statistical accuracy for thermal, 
cold, very cold and probably even for ultracold neutrons. 
In such a case, the asymmetry calculated would be mod- 
ified by an additional short-range interaction as shown in 
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FIG. 6: The bottom two graphs show the ratio of neutron 
flux scattered at an argon atom at an angle of 45 degrees (for- 
ward) to that scattered at an angle of 135 degree (backward), 
as a function of the flnal neutron velocity for the following 
cases: no additional interactions, additional interactions with 
characteristic ranges of 0.1; 1; 10 nm. The thermal motion of 
gas is taken into account. The top two graphs indicate the 
rate of collision with argon atoms (in thermal motion) of neu- 
trons with a fixed initial neutron velocity. The initial neutron 
velocity for the two graphs on the left is equal to 20 m/s; 
the initial neutron velocity for the two graphs on the right is 
equal to 200 m/s. 

fig. ini (the thermal motion of argon atoms is taken into 
account). 

The measurement described above could easily provide 
an accuracy of 10"'^ for the ratio of forward to backward 
scattering probabilities and a corresponding constraint 
for the additional short-range interactions shown in fig.[7l 
The relative drop in sensitivity at a few times 10~^^ m is 
due to the appearance of neutron electron scattering; the 
range of interest for this possible constraint is 10^^^ — 
10-^ m. 



VII. CONCLUSION 

We analysed the constraints for extra short-range in- 
teractions on the basis of the existing data on neutron 
scattering. These constraints are several orders of mag- 
nitude better than those usually cited in the range be- 
tween 1 pm and 5 nm. The reliability of these con- 
straints was supported by the application of several inde- 
pendant methods with comparable accuracy, as well as 
by the use of a major fraction of known neutron scat- 
tering lengths and treatment of the data in a most con- 
servative way. One constraint obtained within the ran- 
dom potential nuclear model was based on the absence of 
an additional linear term in the mass dependance of the 
neutron scattering lengths. It would be difficult to im- 
prove this constraint in either experimental or theoretical 
terms. Another constraint was derived by comparing two 



LIMITS ON EXTRA YUKAWA FORCE 
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FIG. 7: Constraint on extra Yukawa force at 95 % C.L. 
(dashed, dot-dashed and bold lines) obtained in this article. 
The dotted line is an estimation of the sensitivity of the pro- 
posed experiment. Also shown the existing contraints [lll.[l^ 
and the theoretical regions of interest P. Iiol]. 

types of neutron scattering experiments with different 
sensitivities to the extra short-range interactions. These 
are interference experiments measuring forward neutron 
scattering and the Bragg diffraction. The accuracy here 
is limited by the relatively poor precision of the Bragg 
scattering technique. Significant improvements in the 
accuracy of such experiments would be particularly in- 
teresting. Further constraints were estimated using the 
energy-dependence of the neutron scattering lengths at 
heavy nuclei. They are limited by the precision of our 
knowledge of the neutron-electron scattering length. An 
elegant method for further improving such constraints 
would consist in achieving experimental conditions free 
of bne contribution. This is indeed possible, given that 
neutron-electron scattering is essential for fast neutrons 
only. The experiment would consist in scattering very 
cold neutrons at rare noble gases and in measuring pre- 
cisely the differential asymmetry of such scattering as a 
function of the transferred momentum. 
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